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In this work we present a detailed study and derivation of the thermopower and thermoelectric 
coefficient of nano-granular metals at large tunneling conductance between the grains, gr 2> 1. 
An important criterion for the performance of a thermoelectric device is the thermodynamic figure 
of merit which is derived using the kinetic coefficients of granular metals. All results are valid 
at intermediate temperatures, E c 3> T/gr > 5, where 8 is the mean energy level spacing for a 
single grain and E c its charging energy. We show that the electron-electron interaction leads to an 
increase of the thermopower with decreasing grain size and discuss our results in the light of future 
generation thermoelectric materials for low temperature applications. The behavior of the figure of 
merit depending on system parameters like grain size, tunneling conductance, and temperature is 
presented. 

PACS numbers: 73.63.-b, 72.15.Jf, 73.23.Hk 
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I. INTRODUCTION 

Thermoelectric materials with high efficiency is a ma- 
jor research area in condensed matter physics and mate- 
rials science for several decades now. Due to recent ad- 
vances in nano-fabrication, these materials promise next 
generation devices for conversion of thermal energy to 
electrical energy and vice versa. A measure for the per- 
formance or efficiency of a thermoelectric material is the 
dimensionless figure of merit, usually denoted as ZT, 
where T is the temperature. It depends on the ther- 
mopower or Seebeck coefficient, and the electric and ther- 
mal conductivitiesi^^. However, the Wiedemann- 
Franz law has defeated much progress in increasing the 
performance of bulk materials, since it directly relates 
electric and thermal conductivity whereas the figure of 
merit is proportional to the quotient of both. The See- 
beck coefficient of a material measures the magnitude of 
an induced thermoelectric voltage in response to a tem- 
perature difference across that material. If the tempera- 
ture difference AT between the two ends of a material is 
small, then the thermopower, S, of the material is defined 
as S — —Ay/ AT, where AV is the voltage difference 
across a sample. 

In general, thermoelectric devices are used as con- 
verters for either electrical power into heating/cooling 
(Peltier effect) or sources of different temperature into 
electricity (Seebeck effect). These devices are usually 
much simpler, especially without moving parts, than con- 
ventional devices, e.g. two-phase compressors for cooling, 
and therefore more reliable. However, for both effects 
the materials need to have good electrical conductivity 
to minimize ohmic heating and at the same time to be 
bad thermal conductors to avoid thermal equilibration of 
the temperature gradient. Therefore, the aim is to cre- 
ate materials which optimize these parameters together 
with the thermopower. Currently thermoelectric devices 
based on p- and n-type-doped semiconductor junctions 
archive only about 12% of the maximal theoretical effi- 



ciency (as compared to 60% in conventional cooling sys- 
tems )£. 

To be competitive compared with conventional refrig- 
erators, one must develop thermoelectric materials with 
large ZT. The highest figure of merit for bulk ther- 
moelectric materials is about 1, but in order to match 
the efficiency of mechanical systems, ZT ~ 9 is needed^. 
However, for ZT > 2 thermoelectic applications become 
economically competitive^ii^. Although it is possible 
in principle to develop homogeneous materials with that 
large figure of merit, there are no concrete devices on the 
horizon. Especially promising for further improvement 
in efficiency are inhomogeneous/ granular thermoelectric 
materials^ in which one can directly control the system 
parameters. In Ref. [ll[ a figure of merit at 300K of 2.4 
for a layered nanoscale structure and later in Ref. [l2j 
a ZT of 3.2 at about 600K for a bulk material with 
nanoscale inclusions were reported. 

Overall, recent years have seen a remarkable progress 
in the design of granular conductors with controllable 
structure parameters. Granules can be capped with or- 
ganic (ligands) or inorganic molecules which connect and 
regulate the coupling between them. Altering the size 
and shape of granules one can regulate quantum con- 
finement effects. In particular, tuning microscopic pa- 
rameters one can vary the granular materials from be- 
ing relatively good metals to pronounced insulators as a 
function of the strength of electron tunneling couplings 
between conducting grain a 13 ! 14 ' 15 . This makes granular 
conductors a perfect exemplary system for studying ther- 
moelectric and related phenomena. 

All these experimental achievements and technological 
prospects call for a comprehensive theory able to provide 
quantitative description of not only the electric but also 
thermoelectric properties of granular conductors, which 
can in future serve as basis for a clever design of devices 
for a new generation of nano-thermoelectrics. 

Most theoretical progress so far was archived by nu- 
merical solution of phcnomenological model o 16 ' 17 . How- 
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ever, no analytical results obtained from a microscopic 
model for coupled nanodot /grain systems is available up 
to now. Thus, the fundamental question that remains 
open is how thermoelectric coefficient and thermopower 
behave in nanogranular thermoelectric materials. Here, 
we make a step towards answering this question for gran- 
ular metals at intermediate temperatures by generalizing 
our approach^ recently developed for the description of 
electrioi^ and heat transport 2 ^. In particular, we will an- 
swer the question to what extend quantum and confine- 
ment effects in nanostructures are important in changing 
ZT. 

In this paper we investigate the thermopower S, ther- 
moelectric coefficient 77, and the figure of merit ZT of 
granular samples focusing on the case of large tunnel- 
ing conductance between the grains, > 1. With- 
out Coulomb interaction the granular system would be a 
good metal in this limit and our task is to include charg- 
ing effects in the theory. We furthermore restrict our 
considerations to the case of intermediate temperatures 



E c » T/g T > S, 



(1) 



where 5 is the mean level spacing of a single grain and 
E c the charging energy. The left inequality means that 
the temperature is not high enough such that Coulomb 
effects are pronounced. It also allows us to perform all 
calculations up to logarithmic accuracy. The right in- 
equality allows us to consider the electronic motion as co- 
herent within the grains, however this coherence does not 
extend to scales larger than the size of a single graini^. 
In Rcf. 21| we presented a few major results which we 



extend here and describe the derivations in much more 
detail. 

The paper is organized as follows: In Sec. [II] we sum- 
marize our main results and discuss their range of ap- 
plicabilities, in Sec. IIIII we introduce the model, and in 
Sec. IIVI we outline the derivation of thermoelectric co- 
efficient of granular metals in and without the presence 
of interaction which is the main result of this paper. In 
the following section we discuss the behavior of the See- 
beck coefficient and figure of merit (Sec. [Vj as a func- 
tion of sample parameters. Finally, in Sec. IVII we discuss 
our findings and present possible further applications of 
our method. Important details of our calculations are 
presented in several comprehensive Appendixes: in Ap- 
pendix [5] we calculate the thermoelectric coefficient of 
homogeneous disorder metals in the absence of interac- 
tion. In Appendix [B] we derive the heat and electric cur- 
rents of granular metals, and in Appendixes ICl and [Pi wc 
calculate the thermoelectric coefficient of granular metals 
without and in the presence of interaction, respectively. 



thermoelectric coefficient r) of granular metals that in- 
cludes corrections due to Coulomb interaction at tem- 
peratures T > grS, where 5 is the mean level spacing of 
a single grain 
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Here 



(0) _ _ 



(n 2 /3)eg T a 2 - d (T/e F ), 



(2a) 



(2b) 



is the thermoelectric coefficient of granular materials in 
the absence of electron-electron interaction with e being 
the electron charge, a the size of a single grain, d = 
2,3 the dimensionality of a sample, ef being the Fermi 
energy, and E c = e 2 /a is the charging energy. 

The condition for the temperature range of our theory 
ensures that the argument of the logarithm in Eq. (|2"a|) 
is much larger than 1, such that all numerical prefactors 
under the logarithm can be neglected. Furthermore, it 
also defines a critical lower limit for the grain size when 
the charging energy E c becomes of order of the mean 
energy level spacing 5. 

At the temperatures under consideration, the electron 
motion is coherent within the grains, but coherence does 
not extend to scales larger than the size a of a single 
graii>!£. Under these conditions, the electric conductivity 
a and the electric thermal conductivity k are given by the 



expressions- 
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d 



(3a) 
'(3b) 



and = L <7^T, (3c) 



are the electric (including spin) and thermal conductivi- 
ties of granular metals in the absence of Coulomb inter- 
action with Lq = 7r 2 /3e 2 being the Lorentz number. We 
mention that at temperature T > gxS the correction to 
the thermoelectric coefficient, Eq. (|2a|) . has a TlnT de- 
pendence in both d = 2, 3 dimensions which is similar to 
the result for the electric conductivity, Eq. (f3"aj) . having 
a InT dependence in all dimensions as well. 

(ii) Using the above results, we obtain the expression 
for thermopower S of granular metals 



S = (i-l 



Airgrd 



2 , 9tE c 



T 



(4a) 



II. RESULTS AND SUMMARY 

In this section we summarize our results and discuss 
their range of applicabilities. The main results of our 
work are as follows: (i) We derive the expression for the 



where 



S^=-(ir 2 /6)(l/e)(T/s F ), 



(4b) 



is the thermopower of granular metals in the absence of 
Coulomb interaction. 
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(iii) Finally, we find the figure of merit to be: 

Z 7T-2 g T E c 1 f 37, d = 3 

Z(°) 27r. 9T d n T 2TT 2 g T \ ln^, d= 2 ' 

(5a) 

where 

Z(°>T= (7r 2 /12)(T/e F ) 2 , (5b) 

is the bare figure of merit of granular materials and 7 w 
0.355 is a numerical coefficient. In Sec.|V]we present plots 
of Z in dependence of various sample parameters. We 
find, that the influence of granularity is most effective for 
small grain sizes and the presence of Coulomb interaction 
decreases the figure of merit. 

At this point we remark that all results are obtained 
in the absence of phonons which become relevant only at 
higher temperatures. At the end of this paper we will 
briefly discuss their influence. 



III. MODEL 



have an infinite dielectric constant, the effective dielec- 
tric constant of the whole sample can be considerably 
larger than the dielectric constant of its insulating com- 
ponent. Thus the effective single-grain charging energy 
can be much less than the electrostatic energy of a single 
grain in vacuum. 

The last term in Eq. ^ is the tunneling part of the 
Hamiltonian where Uj are the tunnel matrix elements 
between grains i and j which we consider to be random 
Gaussian variables defined by their correlators: 

tfj = tg = const. The dimensionlcss tunneling conduc- 
tance is related to the average matrix elements as 

g T = 2n(t 2 /8 2 ). (8) 

The conductance gx is defined per one spin component, 
such that, for example, the high temperature (Drudc) 
conductivity of a periodic granular sample is = 
2e 2 g T a 2 - d . 



We start our considerations with the introduction of 
our model. We consider a d— dimensional array of metal- 
lic grains with Coulomb interaction between electrons. 
The motion of electrons inside the grains is diffusive, i.e., 
the electron's mean free path £ is smaller than the grain 
size a, and they tunnel from grain to grain. We assume 
that the sample would be a good metal in the absence of 
Coulomb interaction. However, we also assume that the 
tunneling conductance <?t is still smaller than the grain 
conductance go, meaning that the granular structure is 
pronounced and the resistivity is controlled by tunneling 
between grains. 

Each grain is characterized by two energy scales: (i) 
the mean energy level spacing <5, and (ii) the charging 
energy E c = e 2 /a (for a typical grain size of a sa lOnm 
E c is of the order of 2000K) and we assume that the 
condition 8 <C E c is fulfilled. 

The system of coupled metallic grains is described by 
the Hamiltonian H = V\ Hi , where the sum is taken 
over all grains in the system and 



fj t ~\ - 1 Y^ e n ^ n 3 



(6) 

The first term on the right hand side (r. h. s.) of Eq. ^ 
describes the i-th isolated disordered grain, d\ k (&i t k) are 
the creation (annihilation) operators for an electron in 
the state k and = fc 2 /2m— a with /j, being the chemical 
potential. 

The second term describes the charging energy, Cij is 

the capacitance matrix and hi = a\ k a,ik is the number 
operator for electrons in the i-th grain. The Coulomb in- 
teraction is long ranged and its off-diagonal components 
cannot be neglected. Note that, since metallic grains 



IV. THERMOELECTRIC COEFFICIENT 

Having introduced our model in the previous section, 
we come now to the main methodical part of our work, 
the derivation of the thermoelectric coefficient ?/. In gen- 
eral the three kinetic coefficients are: the electric conduc- 
tivity cr, the thermoelectric coefficient 77, and the thermal 
conductivity k. They are related to the Matsubara re- 
sponse functions L^ a ^ with a,/3 € {e, h] 3 i 4 i 23 as 



j(e) = _ 
i(h) _ _ 



(L( ee V(e 2 T))V(eV) - {L^ /{eT 2 )) VT, 
(L( eh V(eT)) V(eV) - (L (Wl) /T 2 ) VT. (9a) 



Here j*- 6 ' (y ) 1S the electric (thermal) current and V 
the electrostatic potential. From Eq. (|9ap one finds that 



(9b) 



a = L^/T, n = L {e ^/T 2 , 
S = -AV/AT = L^/iTL^), 



where the response functions are given by Kubo formulas 



iTd 



1/T 



(10) 

with T T being the time ordering operator for the currents 
with respect to the imaginary time r. Thus, to calculate 
the thermoelectric coefficient rj and thermopower S of 
granular metals one has to know the explicit form of the 
electric j*- 6 -* and thermal currents. 

.(e) 

The electric current through grain i is defined as 



j» e) = $3-4 = edfl *l dt = ie[fk,H] 



>-ifl+5 



(11a) 
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FIG. 1: (color online) Vertices corresponding to the ther- 
mal current operator, Eqs. (|12b[l : vertex (a) corresponds to 



and (b) to j 



•J 



The solid lines denote the propaga- 
tor of electrons, the thick wavy line describes Coulomb in- 
teraction, the tunneling vertices are described by the circles, 
uj n — nT{2n+ 1) and Q. m = 2-kitlT are Fermionic and Bosonic 
Matsubara frequencies respectively (n,m £ Z). 



Straightforward calculations (see Appendix [B| lead to 

(lib) 



Jij 



k,q 



For granular metals the thermal current operator 



•Ah) 



E 



(12a) 




FIG. 2: (color online) Diagrams describing the thermoelectric 
coefficient of granular metals at temperatures T > grS: dia- 
gram (a) corresponds to r]o in Eq. I]2a). Diagrams (b)-(d) de- 
scribe first order corrections to the thermoelectric coefficient 
of granular metals t/ 1 ' in Eq. (|25|l due to electron-electron in- 
teraction. The solid lines denote the propagator of electrons, 
the wavy lines describe effective screened electron-electron 
propagator, and the (red) triangles describe the elastic in- 
teraction of electrons with impurities. The tunneling vertices 
are described by the circles. The sum of the diagrams (b)-(d) 
results in the thermoelectric coefficient correction j/ 1 ' given 
in Eq. J33J). 



can be obtained as follows. The energy content of each 
grain changes as a function of time, such that dHi/dt = 
i[Hi,H]. Energy conservation requires that this energy 
flows to other grains in the system, dHi/dt = Ylj Jij ■ 
Calculating the commutator [Hi ,H], we obtain (for de- 
tails see Appendix IB")) 



Jij 



s(M) 
Jij 



-(h.O) Jh,l) 



*E 



6 + 6, 



L ij u ik a 3Q 



,qk~f - 
ji a jq aik 



E 



c 



.)<>> 



(12b) 
(12c) 

(12d) 



where {A; B} + denotes the anti-commutator. The con- 
tribution v h '°' is the heat current in the absence of 

electron-electron interaction, while the second term 
appears due to Coulomb interaction. Equation (|12b[) im- 
plies that the thermal current operator must be associ- 
ated with two different vertices in diagram representa- 
tion, Fig. [1] We remark that Eq. ^ suggests also a 

finite contribution to jj proportional to t 2 - which in- 
deed exists. However, it vanishes when summed over the 
sample (for details see Appendix [Bj. 

For large tunneling conductance, the Matsubara ther- 
mal current - electric current correlator can be analyzed 
perturbatively in 1 / qt, using the diagrammatic tech- 
nique discussed in Ref. [18| that we briefly outline below. 
The self-energy of the averaged single electron Green's 
function has two contributions: The first one corresponds 
to scattering by impurities inside a single grain while the 



second is due to the process of scattering between the 
grains. The former results only in a small rcnormaliza- 
tion of the relaxation time which depends in general on 
the electron energy to as 



[\ + {d/2-l)u/e F ] 



(13) 



which is a result of the renormalization of the density of 
states at the Fermi surface [see Eqs. (|A8[) and (|A9[) ]. 

In the following we outline the calculation of thermo- 
electric coefficient r\ in the non-interacting case and its 
correction due to Coulomb interaction. A detailed deriva- 
tion of both can be found in Appendices IClandlDl respec- 
tively. 



A. non-interacting case 

First, we consider the thermoelectric coefficient r/ ' 
of granular metals in the absence of interaction. The 
expression for the thermoelectric coefficient in linear re- 
sponse theory is 



,(0) 



d 



)(0) 



(14) 



f2=0 



Here a is the grain size and Q( 0) the correlator of the 
(see Fig. Q}i), and electric current, 



heat current, 



J 



~]( e \ shown in Fig. [2^,. Notice, that there is an impor- 
tant difference between calculations of the thermoelectric 
coefficient rj and thcrmopowcr S and the calculations of 
the electric a and thermal n conductivities in Eqs. pal 
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I3bp . Indeed, to calculate a and k it was sufficient to ap- 
proximate the tunneling matrix element t pq as a constant 
t which is evaluated at the Fermi surface and neglect vari- 
ations of t pq with energy which occur on the scale T/ef- 
However, this approximation is insufficient for calculating 
the thermoelectric coefficient rj and thermopower S since 
the dominant contribution to these quantities vanishes 
due to particle-hole symmetry such that both quantities 
are proportional to the small parameter T/sp- Since it 
is necessary to take into account terms of order of T/e f 
in order to obtain a nonzero result for rj and S, the cor- 
responding expansions must be carried out to this order 
for all quantities which depend on energy. 

For granular metals the important element of the dia- 
gram is the tunneling matrix elements describing the 
coupling between grains i and j. Therefore we derive an 

j 



expression for i.? in the following, assuming that i and 



are nearest neighbor grains and is independent of the 
position in the sample. In order to calculate the energy 
dependence of these elements we assume the tunneling 
barrier between grains to be a delta potential. For the 
one-particle Hamiltonian H = — -^-^ + A<5 (x), with A 
being the strength of the potential, the transmission rate 
for a single particle with energy e p = Ef + £p is 



T — 



2h 2 e 



2h 2 e F 



(15) 



Here Tn 



1X 2 



= T Q {l + T^-i p /e F ) 1 



is the transmission rate at the Fermi 



level and we use the fact that £ p <C £f- Next, we 
consider the case of large barriers, in this regime T p ~ 
To (1 + ^p/sf)- In granular systems we have many chan- 
nels and have to consider tunneling processes with en- 
ergy £j in grain i and with £j in neighboring grain j: 

t 2 (x N (t 2 . + T 2 ^j. Thus, the final expression for tun- 
neling matrix element is 



e F 



(16) 



where t is the constant tunneling matrix clement evalu- 
ated at the Fermi surface. For convenience we use i = 1 
and j = 2 in the following. Using Eq. (JX5J) we obtain 
the following expression for the correlation function in 
Eq. Ull) 



Q 



(0) 



-,2d+2 



(17a) 



f d d Pl d d p 2 



6 +6 



xG(pi,LU n )G(p2,UJn + ftm) • 



where ~n is the unit vector in direction of the current 
a G {e, h}, (lie ■ = l/d is the result of averaging 



over angles, the summation goes over Fermionic Matsub- 
ara frequencies tu n = 2ixT(n + 1/2), and G(p,uj n ) is the 
Matsubara Green's function 



G(p, uj n ) = [vjJ„ - £ p ± i/(2t u )] 



(17b) 



with £ p = e p — £f being the electron energy with respect 
to the Fermi energy [e p = p 2 /(2m)] and the (energy) re- 
laxation time Tuj is defined in Eq. (fl"5|) . To shorten the 
notation, we neglect the momentum argument in the fol- 
lowing and attach the grain index to the Green's function 
G. The two ^-factors under the integration in (|17aj) arise 
from the heat current, Eq. (|12b() . and the energy correc- 
tion of the tunneling clement, Eq. (|16|) . The momentum 
integrals in Eq. (|17a|) arc transformed into energy inte- 
grals taking into account first order corrections in ^/ef 
of the Jacobian [see Eq. (|A7|) of Appendix [5] . 

We first perform the analytical continuation over the 
Fermionic Matsubara frequencies u) n = 2ttT(u + 1/2) in 
Eq. (|17a[) . In order to accomplish that, the analytical 
structure of diagram (a) in Fig. [2] needs to be analyzed, 
which gives rise to three different regions for the Matsub- 



ara summations: I\ = ]— 00; — fi m ] , I2 



-n m ;0[,I 3 



[0;oo[, in which we can determine whether the Green's 
function is retarded, G~(uj), or advanced, G + (lo): 



Si = 



2J G 1 (u> n )G 2 (u n + fi m ) 
n&Ii 



(18a) 



— dio f uj \ _ _ 
4^T tanh V2T J Gl + l ^ G2 (~ IUJ ^ 



S 2 = Gf(bj n )G 2 (cj n + Ct r , 



(18b) 



-duj 



tanh(^) [G^(-tuj)G^(-tu -til) 



S, 



-Gi(-ilj + iQ)G%(-ilj)] 
Y,Gt(uj n )G+(uj n + a m ) 



(18c) 



n£l 3 



' dUJ ■ / Gf, 



tanh — G+ (-iLu)G+ (-VjJ - ity. 



To calculate ry^' 
tive of Si + S 2 -+ 



in Eq. (jT4J) we now consider the deriva- 
S3 with respect to Bosonic frequencies 



Wi\n=o- ^ 0T brevity we omit arguments —ilo of Green's 
functions, leading to 



-diu 



d_ 

\ d 



(Si + S 2 + S 3 ) = (19) 



0=0 



/ ui \ a 
AmT t&nh {2T)d^ 



G-n+n- 
1 1 2 2 



Next, one can perform the integration over variables £1 
and £2 using Eqs. (fT3)) and (|20)) and the residuum theo- 
rem 



/ dfidfa fluCd, ^)GjG+G 2 G+ 



(20a) 



4tt 2 t 2 



2lu + i/t + ^— {2uj + i/t ) 2 
ZEf 
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with 



512 (£1,6) = £1 



■6 + ^ 

q2 2 e F 



(20b) 



As a result we obtain the following expression for the 
derivative of correlation function 



d_ 

on 



)(0) 



0=0 



16« 



-et a 



2 2d+2 



1 



Te, 



(21) 



[2 + (d/2-l)«/ (roej-)] w 



2 , ,2 



— et2a M+2 



3/ 



cosli (w/(2T)) 

2 y2 



Here all contributions of order 1/sp or smaller are ne- 
glected in the final expression. Substituting this result, 
Eq. (|2ip . into Eq. (fl4|) we finally obtain the following ex- 
pression for the non-interacting thermoelectric coefficient 
of granular metals 



V 



(0) _ 



-et a 



2„d+2 



(22) 



One can re-write this expression using the relations: 
vfv d = g T a*- d _ ; ^ = (Sa d )^ ; t 2 = g T S 2 /(2n), 
where T>d is the diffusion constant, gx the tunneling con- 
ductance, and 5 the mean level spacing, giving 



,(o) 



- S -^eg T a 2 - d {T/e F ). 



(23) 



B. correction due to Coulomb interaction 

Now, we consider the correction rf 1 ' to the thermoelec- 
tric coefficient of granular metals due to electron-electron 
interaction 



v = v (o) +77 M 



(24) 



where r]^ is given by Eq. (|22|) . Analogously to ?/ ', j/ 1 -* 
can be obtained from 



r/ 1 ) = t 



0=0 



Q (1) + 



(2) 



Q 



(:!) 



(25) 



where the diagrams Q^ 2 \ contributing to 
are shown in Fig.[2](b, c, d). Detailed calculations of j/ 1 - 1 
are presented in Appendix [Dj However, here we outline 
the main steps of this derivation. These three diagrams 
include the effect of elastic scattering of electron at im- 
purities described by diffusons 



with 



v- 1 = T U (\n i \ + e q S) 



2g F 2d — cos(lf 



(26a) 



(26b) 



where Y^' a stands for summation over all directions and 
orientations |±a~e > ^°' ) j , and the effect of the dynamically 
screened Coulomb potential V(q, f2j) = T>V(q, fli)T> 



V(q,Sli) 

v( q ,n t ) 



2E c (q) 



7-2 00*1+ AE c {q)e q ] [1^1 + e q 6] 

-1 



, (27a) 



1 



2e 



2E c {q) m+e q 5 
where we use the notation 

- ln(qa), d = 1 

EM = - { «/q, d = 2 

2n/q 2 1 d = 3. 



(27b) 



Each diagram in Fig.[5]has also two types of renormalized 
interaction vertices: (i) the inter-grain vertex 



a d d~g 2E c (q)J2' a cos(~q ■ ~a) 
(2n) d T* [|n,.|+4£7 c (g)e g ] + e q S] 

and (ii) the intra-grain vertex 

a d d~q 2E c (q)2d 



(27r) d T-2 [\Sli\ + 4E c (q)e q ] [\SU\ + e q S] 



(29) 



Explicitly, the contribution Q^> in Eq. (|25|) [diagram (b) 
in Fig. [2] is given by 



Q 



(i) 



2d 



-et 2 T 2 a 2d+2 



(^ 0, ) ; 



(30a) 



where the function g 12 is defined in Eq. (|20b|) and we use 
the notation 



(30b) 



with Sj = ± denote the analytic structure of the Green's 
functions implying restrictions on the frequency summa- 
tion. 

For the contribution in Eq. (f25|) [diagram (c) in 
Fig. [2] we have the following expression 



Q(2) = ^-L et 2 T 2 a 2d+2 (^ 0) J (31a) 
xE Jd^ 29l2 Ft S2S3) ^ 2) (^), 
where we use the notation 

(31b) 

The diagram Q^ 3 - 1 , shown in Fig. [2] (d), describes the 
contribution of the correlation function with the interac- 
tion part of the heat current operator, J^' 1 ^ [second term 
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FIG. 3: (color online) Plots of the dimensionless figure of 
merit Z/Z^- ' vs. grain size a (in nm) - for different values 
of dimensionless tunneling conductance gr (see legend): the 
upper panel is for the three-dimensional (3D) case and the 
lower for the two-dimensional (2D). All curves are plotted for 
T = 100K. At this temperature, the dimensionless bare figure 
of merit for granular metals is Z^T ~ f0~ 4 . 
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FIG. 4: (color online) Plots of the dimensionless figure of 
merit Z/Z^ ' vs. tunneling conductance gr - for different 
values of grain sizes a (see legend): the upper panel is for 
the three-dimensional (3D) case and the lower for the two- 
dimensional (2D). All curves are plotted for T = 100K. 



in the right hand side of Eq. (|12bj) ] , and has therefore a 
different structure in comparison with contributions 
and Q^: 



(3) 



Q 



where 

(sis 2 s 3 ) 



2d et l a 



2rp2 2d+2 I ,,(0) 



F. 



' E / < '^'i^ * 3 («<>«) 



f]?, 



2 1 



d 
2e~f~ 



& + 6 



(32b) 



The main contribution to r/ 1 ) from diagram is of 
the order of (T/ef) 2 , whereas and have 1/ef 
contributions. Therefore we will not consider diagram 
any further, but keep contributions of order T/ef 

only. 

Thus, the first order interaction corrections to the ther- 
moelectric coefficient are only generated by diagrams (b) 
and (c) in Fig. H Substituting Eqs. (gOg) and (|3Ta|) 
into Eq. (|25p after summation over Fcrmionic, w n , and 
Bosonic, f2j, frequencies and analytical continuation we 
obtained (see Appendix [D] for details) 



n 



(i) 



2ngT \2tt 



) / d d qln 


'2D c (q)e q " 




T 



(33) 



where the q-integration goes over the cf-dimensional 
sphere with radius it /a. Integrating over q in Eq. (|33|) 
we obtain the following expressions, neglecting all con- 
stants under the logarithm, in two (2D) and three (3D) 



dimensions 



V (0) ,„ E c g T (1 



In- 

?t T 



%D 



■In 



T : 



(34) 



(32a) which lead to Eq. (2a] 



V. THERMOPOWER AND FIGURE OF MERIT 

Now we have the expressions for all three kinetic coef- 
ficients (T, K, and rj for granular metals to order T/ef- 
Based on these, we can derive other thermodynamic 
quantities - in particular we discuss the thermopower 
(Seebeck coefficient) and figure of merit in this section. 
Both quantities are relevant parameters for thermocou- 
ples, the thermopower is a measure the change of voltage 
due to a temperature gradient and the figure of merit a 
measure for the performance of the device. 

The thermopower is related to the kinetic coefficients 
as (see Eq. (|9bl ) 



S = r\ja . 



(35a) 



Again, we only consider terms up to order T/ef and 
obtain the expression 



77W /77W 



(35b) 



which results in Eq. (|4a[) . 

The dimensionless figure of merit is related to the ki- 
netic coefficients as 



ZT = T7] 2 /(<tk) = S 2 aT/n, 



(36a) 
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graphs) and gr — 2 (lower two graphs) and grain size 
a — 5, fOnm and d = 2,3. Legends are next to the corre- 
sponding curves. 



giving 



(s(°))V°)r / ;cW qW 2 gw\ 



Zr= ^ " l 1 -^) + ^T + 2 ^))' (36b) 



resulting in Eq. (|5a|) . which has lowest order (T/ef) 2 
[Eq.EE]. In Eq. |36b|) qW, and S^ 1 ) are corrections 
to the thermal conductivity, electrical conductivity, and 
the Seebeck coefficient due to Coulomb interaction, re- 
spectively. The numerical coefficient two in front of the 
last term reflects the fact that the Seebeck coefficient ap- 
pears squared in the definition of ZT. Using Eqs. (pEaj) 
and (|36bp one can see that the second term of the right- 
hand-side of Eq. (|5a| originates due to correction to the 
Seebeck coefficient, Eq. ([4"af. The origin of the third term 
in Eq. (|5a|) is due to correction to the thermal conductiv- 
ity, K, last term on the right hand side of Eq. (|3b|) . Wc 
mention that the second term on the right hand side of 
Eq. (|3b[) cancels with the correction to the electrical con- 
ductivity, Eq. (|3"a|) . after the substitution into Eq. (|36bl) . 

In Fig. [3] the dependence of Z on the grain size a and 
in Fig. 3] on the tunneling conductance gr for two- and 
three-dimensional samples are shown. Fig. [5] shows the 
temperature dependence of the figure of merit. These 
plots show that the correction term to ZT is most effec- 
tive for small grain at not very high tunneling conduc- 
tance and at low temperatures. 



VI. DISCUSSIONS 

In the presence of interaction effects and not very low 
temperatures T > grS, granular metals behave differ- 
ently from homogeneous disorder metals. However, in the 
absence of interactions the result for rf ^ below Eq. (|!?a|) 
[or Eq. IC15j coincides with the thermoelectric coefficient 
of homogeneous disordered metals (see Appendix [XJ, 



with pp being the Fermi momentum. One can expect 
that at low temperatures, T < grS, even in the presence 
of Coulomb interaction the behavior of thermoelectric co- 
efficient and thermopower of granular metals is similar to 
the behavior of r\hom and Shorn , however this temperature 
range is beyond the scope of the present paper. Our re- 
sults for thermopower (|4ap and figure of merit (|5ap show 
that the influence of Coulomb interaction is most effective 
for small grains. The thermopower S 2 decreases with the 
grain size which is a result of the delicate competition of 
the corrections of thermoelectric coefficient (p?aj) and the 
electric conductivity (|3a[) . In particular, if the numerical 
prefactor of the correction to rj would be slightly smaller, 
the sign of the correction to S would change. 

Above we only considered the electron contribution to 
the figure of merit. At higher temperatures T > T* , 
where 



\j9TC 2 ph /l ph a, 



(38) 



is a characteristic temperature with l pn and c p h being 
the phonon scattering length and phonon velocity respec- 
tively^, phonons will provide an independent, additional 
contribution to thermal transport, 



^ph 



T 3 / , I r 2 



(39) 



However, the phonon contribution can be neglected for 
temperatures 



g T S<T < T* 



(40) 



(o) 

^hiom 



-(2/9)ep F (r T), 



(37) 



A detailed study of the influence of phonons at high tem- 
peratures, including room temperature, will be subject of 
a forthcoming work. 

So far, we ignored the fact that electron-electron in- 
teractions also renormalizc the chemical potential fi. In 
general, this renormalization may affect the kinetic coef- 
ficients: the thermal current vertex, Fig[TJ as well as the 
electron Green's functions depend on fi. In particular one 
needs to replace V(eV) -> V(eU+^) in Eq. ([9"a|), To first 
order in the interactions, the renormalization of /i only 
leads to corrections to diagram (a) in Fig. [2j As it can 
be easily shown, for this diagram the renormalization of 
the heat and electric current vertices is exactly canceled 
by the renormalization of the two electron propagators. 
Therefore, the renormalization of the chemical potential 
does not affect our results in the leading order. 

Finally, we remark that the bare figure of merit Z^T 
for granular metals at gT > 1 and I00K is of the order 
of only I CU 4 . Therefore these materials are not suitable 
for solid-state refrigerators, but should be replaced by 
granular semiconductors with gx < 1. However, the case 
of granular metals is still relevant for low temperature 
applications in, e.g., thermocouples. Therefore we con- 
clude this paragraph discussing the dimensionless figure 
of merit ZT of granular materials at weak coupling be- 
tween the grains, gx <C 1. In this regime the electronic 
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where the diagrammatic representation of the correlator 
A®' is shown in Fig. [5] with electric and heat current 
vertcxes 



Ji 



e 
m 



Jh 



1 

m 



^P^papap, (A2) 



FIG. 6: (color online) Lowest order diagram for the heat- 
electric current correlator for homogeneous disordered metals. 
The external Bosonic frequency is denoted by Q (wavy lines) 
and the internal Fermionic frequency by uj (straight lines). 
The electric and heat current vertexes are j e and jh, respec- 
tively. 



contribution to thermal conductivity ««■ of granular met- 
als was recently investigated in Ref. \24 . where it was 
shown that 



g 2 T T z /E 2 c 



(41) 



In this regime the electric conductivity of granular metals 
obeys the la w 18,22 



cr ~ g T exp(-E c /T). 



(42) 



However, an expression for the thermoelectric coefficient 
in this region is not available yet, but recently it has 
been proposed, based on experiment, that nanostruc- 
tured thermoelectric materials in the low coupling region 
(AgPb m SbTe 2+rn , Bi 2 Te 3 /Sb 2 Te 3 , or CoSb 3 )^^^ 
can have higher figures of merit than their bulk counter- 
parts. 

In conclusion, we have investigated the thermoelectric 
coefficient and thermopower of granular nanomaterials 
in the limit of large tunneling conductance between the 
grains and temperatures T > gxS. We have shown to 
what extend quantum and confinement effects in granu- 
lar metals are important in changing ZT depending on 
system parameters. We also presented the details of our 
calculations. 
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APPENDIX A: THERMOELECTRIC 
COEFFICIENT OF HOMOGENEOUS 
DISORDERED METALS IN THE ABSENCE OF 
INTERACTION 

In order to demonstrate important steps of our cal- 
culations, we present a derivation of the thermoelectric 
coefficient for homogeneous disordered metals in the ab- 
sence of interaction in this appendix. In linear response 
theory the thermoelectric coefficient can be written as 



,«>) 







L d Tdn 



A (0) 



(Al) 



The sums over momentum arc 

d Jd d p by 



, 27W 



respectively. 

transformed into integrals ^2 

p 

which we can transform the sum over the mo- 
mentum product to J d d p e J d d ph p e • P h = 

j d d Pe [d d Ph n? e | 2 07? e • n h ) 6 o? e - y h ) (^) d , 

where n a is the unit vector in direction of the current 
a G {e, h}. Averaging over angles gives ( n e ■ ~n h) = 1/d, 
and therefore: 



a m z J (2tt 



d d p .-.a 



x£ p G(p + k,u n + {l m )G(p,uj n ), (A3) 

where s is the spin degeneracy factor, fl m an external 
bosonic Matsubara frequency, and G(p,ui n ) the momen- 
tum dependent Matsubara Green's function. In the fol- 
lowing we only consider the case of zero external momen- 
tum, k = 0. For the advanced (in C + ) and retarded (in 
C~) Green's function, we use 



G ± (jp, u>„) = [iu n -£ p ± i/{2t^ 



(A4) 



where £ p = e p — Ef is the electron energy with respect 
to the Fermi energy [e p = p 2 /(2m)] and the (energy) 
relaxation time r w depends on the (real) frequency ui. 

The momentum integral in Eq. (|A3[) is transformed 
into an energy integral as follows 



d d p 



(27T 



3/(171) = 



m -2 



(2ir) d 



/(v / 2me)e2 1 de, 
(A5) 



where fid is the value of the angular integral (0:1,2,3 = 

{2, 27T, 47r}). Since the Green's function depends on £, 

we need to rewrite this integral using the de = d£ and 

/2-1 



£ d/2~l = £ d/2-l 



F + 

e d / 2 - l de = 4 /2 - 

d/2-l 



/OO 
[l + (d/2-lK/e F ]^- 
-00 



Combining Eqs. (|A5j) and (|A6|) we obtain 



d d p 



(2n 



-J(\P\) 



(0) 



with v 



(0) 



l n(d) 



y d — (27r) d F 

states (DOS) at the Fermi surface 



(A7) 



d l 2 e d J 2 1 being the density of 



xf(^/2m(Z-s F ))dZ 
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The DOS (and r w ) depends on u> via 



(A8) 



,(0) 



1 



e F 



1 + x 



jx) d l 



2-1 



1 



-dx , 



with x = 2r w £i?(l + lu/ef) and the symbol stands for 
imaginary part. However, the term x/x is neglected in 



the integral, hence [L 



l) dx — 7r] 



(A9) 



For conveniens we drop the momentum argument of G 
in the following, as well as the — iuj argument. 
Using all the above equations we finally obtain 



(o) 2sie (o) 9 



d m 



an 



i2=0 



(A10) 



x[l + {d/2-l)Z/e F ]Z(Z-e F ) 
1 

^ 

i((jj„ + Cl m ) — £ + 1 sgn(w„ + fl m )/ (2r) 
1 



iu) n - £ + isgn(w„)/(2r) ' 

Here the derivative by the real external frequency, Jl, 
requires the analytic continuation of the Matsubara ex- 
pression. For the following calculation of the ui n sum and 
energy integral, it is convenient to introduce the notation 

9i = [l + {d/2-l)t/e F ]t(t-e F ) (All) 

= -& F + (2 - d/2)e + o [e/s F ] . 

Next, we split the sum over the Fermionic Matsubara 
frequencies u> n = 2TrT(n + 1/2) into three intervals I\ = 
]— oo ; — f2 m ] , I2 = ]— ^m! 0[ , ^3 = [0; 00 [ such that we can 
specify the analytical structure of the Green's function, 
i.e. 



2s le 



(0) 



8 



0=0 



d m d dn 

neli 
nel-i 

- J2G + (io n + {l m )G+ 



Now, we perform the analytical continuation of the sums 
to real frequencies and at the same time of the external 
frequency: 



Si = ^2 G (ui n + O m )G (w„) 

neli 

luJ : t£inh(-^-)G~(-uj)G~(-tLJ + ifi) , 



(A12) 



4mT 



'2T 



S 2 = ^G+K + a n )G-K) = J _^tanh(^)(A13) 
x [G + {-iuj - in)G-(-iu) - G + (-ilu)G- (-leu + ifl)] , 



S3 = ^G+K + fi m )G + (w„) (A14) 
^tanh^G+H^-z^G+t-^). 



Next, we rewrite the il-derivative of the integrant in 
terms of ^-derivatives, (arguments — ius are omitted): 



d_ 

an 



0=0 



ouj 



(A15) 



which is valid for both, advanced and retarded functions, 
and therefore we can simplify 
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d_ 

an 



'-G+i-tu - tfl)CT + G+G-(-ilu + ifi)] = (G + ) 2 G~ + G+ (G~) 2 = - A (G+G~) , 



(A16) 



_9_ 



[-G-G-{-iu + iil)] = 



d_ 



[G + {-ilu - iQ)G + ] = 



fi=0 



Using Eqs. (|A15|) and (|A16|) we can write the following 
9 



[5i + S 2 + S 3 ] = 
0=0 



8ttiT 



We notice that at this point the w-integration should not 
be done by parts, since the boundary term is important. 
We can now do the ^-integration using the fact that G + — 
G~ = -(i/t)G+G- 



5( W ) 



d Z 3? 



iG+G~ 



(A18) 



(r[ W -€ + »/(2r)][ W -€-«/(2r)])' ! 



-2tt 



rw 2 (4-rf) + 



4r 



Here, the higher order terms in are neglected. In 
Eq. (|Al8j) we need to keep the terms proportional to oj 2 
only therefore after w-expansion of t = t w we obtain: 

E(cj) = -2tt [t lu 2 (A -d) + Lu 2 T (d - 2)] = -Aitt u} 2 . 

(A19) 

And. finally 



s lev. 
d m 



(0) 
d 



duj , a; . d . 
tanh — )— — s (w) 



s e jo) T 



2d to 



cosli (w/2T) 



2it2s e (0), ^ 



(A20) 



3d to 

To perform the last integration we used the integral 
/ cosh^Or) = tt 2 /6. In d = 3 the density of states has 
the form v^l z = ^rfl 1 , leading to 



(o) 
Wd 



9 



epF(r T) . 



(A21) 



APPENDIX B: HEAT AND ELECTRIC 
CURRENT OPERATORS OF GRANULAR 
METALS 

In this Appendix we derive an expression for the heat 
and electric current operators of granular metals in the 



presence of Coulomb interaction. The Hamiltonian for 
the granular system is H = Y^ej where 



k a i,k 



k a i,<l + l ji a j,q a i,k 



(Bl) 



Here we introduce the notation n, = Yla] jfii.k for the 

k 

number of electron within a grain i. The creation (annihi- 
lation) operators a| k {fli^k) satisfy the anti-commutation 

relations {at fc ;Oj l5 } = 8 l] 8 kq and ja$;o^}_ 

The electric current through grain i is 



0. 



drii 
l ~dT 



o ^ — ^ 



(B2) 



And the heat current through grain i is 



dei 
i—— 
dt 



efjH 



(B3) 



Accordingly the total current operators are given by 



i,3 



(B4) 



where are the unit directions of the current (electric 
field, E, gradient) or heat (temperature, T, gradient) 
flows. 



1. electric current operator 

First, we calculate the electric current. Since the num- 
ber operator of electrons n, commutes with the first two 



terms of H in Eq. (|B2[). i.e. 
need to calculate 



nil H 



;W<*> 



, we only 
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fc' 



k,i' ,j,k' ,q 



(B5) 



r 



As a result we get the following expression for the electric 
current operator 



fii = le 



£ ifiA 

k,q 



k^Jtl "ji a j,q L 



(B6) 



2. heat current operator 



The heat current operator can be conveniently written 
as a sum of two contributions 



Jij Jij 



where, the non-interacting part of the heat current 
originates from the sum of the commutators cj et ^+ c\ te ^ 



(B9) 



o j , , , , , c i-i, i and the interacting part fa from the sum of cj ct - 

becond we turn to the heat current, for which we need , ° ^ J w 1 



to calculate nine commutators (see Eq. (|B3[) 1 



C 



(B7) 



where a € |e,c, ij and operator was defined m „, „ . , , . , ,. ... . 

„ jtvttt 1 We first calculate the non-mteractmg part of the heat 

Eq. (|Bip . ^f/i oi 

Since operators commute with each other, four of current operator, >. That is, we consider the sum of 

the commutators in Eq. (|B7[) vanish: commutators Cf + C\ \ . A straightforward calculation 



a. non-interacting part of heat current operator 



(B8) leads to 
I 



cS et) 



rite) 



1 

2 

i' .j,k' ,q.k 
j,k,q 

\ £ ^ 

j,k,q,i' ,k' 



— tjj a j, q a i,k 



'ij u i,k a J,<l ^ l ji u j,q L 



a], k ,ai>,k' —a\, k ,ai' t k' 



/. . A- 



j,k,q 



(BIO) 



(Bll) 



r 



Using Eqs. (|B10|) and (|B1 1[) we obtain for the sum 



As a result the non-interacting part of the heat current 
operator of granular metals has the form 



;/■'/ 



(Zk+Qtpl q a i>k ). (B12) 



^' 0) = 



£ 



& + e 9 



t kq tf a - t qk a^ a u 

L ij u i,k u 3,<l ''ji a j,q l ' k 



(B13) 
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b. interacting part of heat current operator 



k,q 



Using the calculation for 



To obtain the expression for the interacting part of Cj [Eq. (|B5[)] we get 
the heat current operator we need to consider 



the commutators of the form 



■At) ~ 



where ef) 



At) ~ 



E Vi'j { s ij s *«M> ~ Su'Skk'alkHgj + tf v \5 w 5 kk ,a)^i, k - Mfc^S^)) (B14) 



k : k' ,q 



k.q 



\Ae) 



For the following steps of the calculation of C\ 1 and C\ \ we need the commutator 



k,q 

- (S m j - 5 m i) + ^oi^.fe) = 2 (<*mi ~ <W ? 



At) 
ij ' 



(B15) 



To calculate the interacting part of the heat current operator we need the following commutators 



^nitc) 
,2 i 



2Y,(*C^%$], = ; E ^ (W - + fe' - <W (B16) 



3,i',j' 



I Jrrn Jim J "j ' "t IJmj Jjtn 



2^ 

i 

J, ™ 

^Et^'*-^^^) = 4 E ^ - + - ( B17 ) 
-£ 

3 

~gE [(^mj ~ ^mi) 



ji v ji 
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Using the symmetry relations: = —J^ and C m ] = C^, and 



- .He) 



2ie (S ml - Sm^efj [Eq.|B15J we obtain 



le ' 
1 



3,m 



^{e)~ ~ ^(e) 
3mi n j ' n iJmj 



^mj y ™i Jij j + ^mi \Jij ' "™ J 



; E 



3« i» »j J» 



= 



(B18) 



As a result for the interaction part of the heat current operator of granular metals we obtain 



Jij 



4^ 



(B19) 



So far we have omitted the last commutator in Eq. (|B7[) . which would be an additional contribution to the 
non-interacting part of the heat current operator, 5# m Eq. (|B13j) . However, if this term is summed over i it 
vanishes and docs not contribute to ■ However, for completeness we present the calculation of Q here as well. 

We define the commutator C) as 



j,i',j' 



A-B, 



(B20) 



and simplify first the expression for operator 4B. 



i'j'ij 



a i',k' a r,q' + tf'i' a j- A - a V.v) [ t ij a i,k a 3,9 + l ji 



k,q,k' ,q' 



k,q,k' : q' 

I fk'q'.gk 

U i'j' L ji 



.q'k'.kq 
b ij 



^j'j^q'q a i',k' a i,k ~ ^i'i^k'k a j,q a j',t 



.q'k'.qk 



&i'i&k'kO>j' , q iO>j,q ~ ^j'j^q'q a \.k a i',k' 

Si'j6h' q a! j , q ,ai t k - SfiSq'kaj !q ai',k 



(B21) 



AA. 
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And the final step is: 



--E 

j,k,q 



" J» U j,q"" 1 >k' 



"m.p a j,q • l mi l ij a m,p a j,q 



- E (##4,^ 

j,m,k 7 q 7 p 

" mj ji a m,p a i,k "+~ l mj l ji a m,p a i,k ~ l i m l ji a j,q a "m,p ~ l i m l ji a j,q a 
2 | l rni l ij a m,p a j,q ~ l jm l ij a i,k am ,P + l mj l ji a m,p a i,k 



l 3'rji a j',q' a hk 



mp t kq~\ - 



,qk^ 
l im l ji a j,q a m,p 



j,m,k,q,p 



\ E [ 

j,m,k,q,p 



,qp ,kq^ ^ ,pq ±qk<rfi 



t q Lt': q 



j,m,k,q,p 



®i k^^-V ®tn,p®i>k 



tiL which is 



The underlined terms cancel since ijL 
also used in the last step. If this expression is summed 
over i, we can exchange indices m and i (and A: and g) in 
the second summand. As a result we obtain 



£5 



0. 



(B23) 



Q<°> 




i.e. there is no additional contribution to the full heat 
current operator j h in Eq. (|B4|) from this term. 



3. summary: electric and heat current operators 



To summarize this appendix we explicitly write the ex- 
pressions for electric jjj and heat j\j current operators 



-(e) 
7- 



Jij 



Jij 



"E (*« ^,fc« 



J, 9 



(B24) 



(B25) 
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''ij a ik Uj 31 
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APPENDIX C: THERMOELECTRIC 
COEFFICIENT OF GRANULAR METALS IN 
THE ABSENCE OF INTERACTION 

In this appendix we consider the thermoelectric coeffi- 
cient t/ ' of granular metals in the absence of interaction 



FIG. 7: (color online) Lowest order diagram for the heat- 
electric current correlator of granular metals. The external 
Bosonic frequency is denoted by Q (wavy lines) and the inter- 
nal Fermionic frequency by ui. The electric and heat current 
vertexes are j e and jh, respectively. The ± denote the pos- 
sible analytical structure of the Green's functions (straight 
lines). 



in analogy to Appendix [AJ The expression for the ther- 
moelectric coefficient in linear response theory is 



,(o) 







)(0) 



(CI) 



Here a is the grain size and the correlator of the heat 
current. (see also Fig.[T^,), and electric current, ~j e , 
shown in Fig. [7] 

For granular metals the important element of the di- 
agram is the tunneling matrix elements t*j? describing 
the coupling between grains i and j. Therefore we de- 
rive an expression for t^j in the following, assuming that 

i and j are nearest neighbor grains and is indepen- 
dent of the position in the sample. In order to calculate 
the energy dependence of these elements we assume the 
tunneling barrier between grains to be a delta potential. 
For the one-particle Hamiltonian H = — + X5 (x) 

the transmission rate for a single particle with energy 
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£ P = £f + 6 is 
mX 2 



1 



2h 2 e 



mX 



= T„ 1 



2h 2 e F 

2h 2 e F 
m\ 2 



(I-W^f)J (C2) 
i 

- 6/ £ f 



Here T 1 



is the transmission rate at 



£f and we use the fact that 6 <C £f- Next, we con- 
sider the case of large barriers, in this regime T p ~ 



To (1 + ^p/sf). In granular systems we have many chan- 
nels and have to consider tunneling processes with en- 
ergy £i in grain i = 1 and with £2 hi grain j = 2: 
t 2 oc N (T 2 ± + T 2 2 ). So, we obtain 



n6,6) 1+ 



6 + 6 



(C3) 



Therefore we have the following expression for correlation 
function in Eq. (|C1[) 



Q 



(0) 



2d 



et 2 Tfl 2d + 2 (JO) 



(o)\sr 2^+2 /" rf V f d d P2 ( 6+6 \ A , 6 + 6 A ^ ^ , x 

E / ^1^2 ff(6,6)G(fc,w n )G(6,w n + fi n 0, (C4) 



where ~n is the unit vector in direction of the current 
a G {e, /j}, (^e "* • 'n'h^') = 1/d is the result of averaging 
over angles, the summation goes over Fcrmionic Matsub- 
ara frequencies w n = 2nT(n + 1/2), and G{p,ui n ) is the 
Green's function defined in Eq. (|A4|) of Appendix [Al with 
momenta/energies Pi/6 of grain i. To shorter the nota- 
tion in the following we neglect the momentum argument 
and attach the grain index to G. In Eq. (|C4[) we intro- 
duce the notation 



<?(6,6) = (6+6) 1 + 



6+6 



X 1 + 



d --l 

2 



^111 



4-1 
2 



Si? 



(C5) 



r 



The factors in this order arise from: the heat current 
vertex Eq. (|12cp . the energy correction to tunneling ele- 
ments Eq. (|16p. and corrections to the DOS due to finite 
Fermi energy Eq. (|A7|) . In the linear order in £/ £ f we 
obtain 



3(6,6 = 6 + 6 + o 

2 Ef 



(C6) 



We first perform the analytical continuation in 
Eq. (|C4[) (for convenience the grain index is written as 
index to the Green's functions) 



dtidt* 3(6,6) 



G 1 K)G 2 K + n m )+^2 G+{uj n )G 2 ( 

)+ (Wn)G 2 + (w„ + fi m ) 

ri£/i ra£/2 n£/3 



Si 



S 2 



After analytical continuation we obtain 

5l = ~/ 4^T tanh (^) G rH w +^) G 2(-^), 



s 2 = - 
S3 = - 



_tanh(|f) [Gr(-tw)G+(-K«;-*fi)-G; (-10; + tfi)G+(-tw)] 



4ttiT 



tanh 



(2^) Gt(^u)Gt(-iu-in). 



(C7) 



(C8) 
(C9) 
(CIO) 



Now wc consider the derivative of the sum of (Si + S2 + S3) with respect to Bosonic frequency ^ \q—q- ^ or brevity 
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wc omit arguments — ilj of Green's functions 

{Si + s 2 + s 3 ) = [ - dl 

f2=0 



d_ 

on 



^ tanh ( — — ) 
AmT \2T) 



^gAg 2 -g x (£g+ 



du> ( lu \ 

— tanh — — 

AniT \2TJ 



exchange indices 



dtU j LU 

— tanh — - . 

AniT \2TJ 



c 2 - g£) (^i - ^1 ) 



1 

—2^1 1-^2 t - z 2 



(Cll) 



1 2[{Gi-G+)^{G--G+) + {G--Gt)M G 2- G t)] 

- I £t (w) I V* - «0 W -«)]"/ ^ -* (w) I 

Now one can perform the integration over variables £i and £2 using Eqs. ()C7|) . (|C11|) . and the residuum theorem 

/ d^ 2 g(^,^)G^G+G 2 G+ (C12) 



4i +42 + 5 — — 



(w - fi - i/(2r)) (w - ix + i/(2t)) (w - 6 - i/(2r)) (w - f 2 + i/(2r)) 



2ttc / d£ 2 



Ci,o=w+«/(2t) 

o; + »/(2r) + + ^ (tj + e/(2t) + &f 
i-t/r) (u - 6 - i/(2r)) (w - 6 + */(2t)) 

" V ' 

f2,0=W+j/(2T) 



47r 2 r 2 



2lo + i/t+^— (2uj + i/rf 

IEf 



where the 4^0-values below the braces denote the poles in the complex plane used to perform the integration. As a 
result we obtain the following expression for the derivative of correlation function 



d_ 

on 



lTd et ° a 



-™etW d+2 (vf) (2T)^-^ / duo 



4?r 2 2lo + i/t + — (2lu + i/t) 2 
2ef 

^-{2uj + i/T a [{d/2-l)u:/e F ]f 

ZEf 



(C13) 



Aid 



= -—et n a 
Aid 



2 2d+2 / .(0) 



w (2TY— — . 
d 1 v ' 2e F 6 



\2T J dto _ 
1 d r d J2 + (d/2-l)i/ (t o£f )) 2 ^ 2 
2e F J ^ cosh 2 (w/(2T)) 

4d 7T 2 



In the second line in Eq. (|C13j) the derivative is taken 
into account (removing the boundary terms of the par- 
tial integration) and the contributions of order 1/sp or 
smaller are neglected in the last line. 

Finally, we obtain the following expression for non- 
interacting thermoelectric coefficient of granular metals 

,(0) = _f|! e ^ +2 ^(0)) 2 ^. (C14) 

One can re-write this expression using the relations: 
vfv d = ga 2 ~ d ; vf = (S^r 1 ■ tl = 9 6 2 /(27t), where 
T>d is the diffusion constant, g the tunneling conductance, 



and S the mean level spacing, giving 

V {0) = -^-egW- d {T/e F ). (CIS) 



APPENDIX D: THERMOELECTRIC 
COEFFICIENT OF GRANULAR METALS WITH 
INTERACTION 

In this appendix we consider the correction rf 1 ' to 
the thermoelectric coefficient of granular metals due to 
electron-electron interaction, i.e., 

r? = r? (0) +7? (i) ) ( D1 ) 



18 




FIG. 8: (color online) Diagram describing the correction to 
the thermoelectric coefficient due to electron-electron interac- 
tion corresponding to the term Q' 1 ' in Eq. (|D2[) . The external 
Bosonic frequency is denoted by Q. (wavy lines) and the inter- 
nal Fermionic frequency by u) (straight lines) . The electric and 
heat current vertexes are j e and jh (without Coulomb contri- 
bution), respectively. The (red) triangles denote the diffusons 
T> and the thick wavy line the screen Coulomb interaction. 



where j/ - 1 was calculated in Appendix [Cl Eq. (|C14|) . The 
structure of the diagrams 

Q (2) , Q {3} contributing to 
r/^ are shown in Fig. [5] and we can write 



a d Td£l 



g (1) + q w + q 



(2) 



)(3) 



(D2) 



These diagrams include the effect of elastic scattering 
of electron at impurities described by diffusons, 2? _1 = 
t uj (l^il + £(j<5), and the effect of the dynamically screened 
Coulomb potential V(q, = T>V{q, Vt^V: 



V{q,SU) 
V(q,ni) 



2E c {q) 



Tl[\ni\+AE c {q)e q }[\^i\+e q S\ 
l 



, (D3) 



1 



2e„ 



2E c {q) p i \+e q S / 
The renormalized interaction vertices arc: (i) inter-grain 

W 2T- 2 E c ( q )J2' a cos(li -1?) 



(2w) d [\n t \+4E c (q)e q }[\n t \+e q 5y 



and (ii) intra-grain 

$£ 2) (a) 



(D5) 



d d dq 



4dr~ 2 E c (q) 



{2Tr) d [\SU\ + 4E c (q)e q ] [\ni\+e q Sy 



with e q = 2gx (2d — ~Y^ a cos(~if • ~a)) where ^2' a stands 
for summation over all directions and orientations 
|±a"ej°'|, and 



e 2 e 2 ( -ln(ga), 
2C(q) a d 1 27r/g2j 



d=l 

d = 2 (D6) 
d = 3. 



Explicitly, the contribution Q^ 1 ' in Eq. (|D2|) is given by 



Q 




FIG. 9: (color online) Diagram describing correction to the 
thermoelectric coefficient due to electron-electron interaction 
corresponding to term in Eq. (|D2[) . All notations are the 
same as in Fig. [8] 



-,(D 



s 

Yd 



et 2 T 2 a 2d+2 ^ (D?) 



Q (1,±± 



with gi2 = 5(^1,^2) defined in Eq. 
Here we introduce the function 



C6|) of Appendix O 



(S1S2S3S4) 



GJ 1 (w„ + ft m )G* 2 (w„ + ro + 0;) 



x G s 2 3 (w n )G| 4 (w„ + f2i), 



(D8) 



with Si denote the analytic structure of the Green's func- 
tions implying restrictions on the frequency summation 
- in principle there are 16 different combinations of the 
Si's, see Fig. [8] However, only the two diagrams <3( 1-±± ) 
and Q( 1,=F=F ) contribute to the correction rf 1 ^ . The other 
analytical structures do not have either valid frequency 
domains or poles are located in only one half plane of C. 
For the contribution in Eq. (|D2|) . Fig. [H we have 
the following expression 



)(2) = _±- et lT 2 a 2d+2 



,(0) 



E 



Q 



(2,±+) 



Q 



(2,T- 



(2,±-) 



where 



(D10) 

Due to symmetry all three contributing diagrams in the 
right hand side of Eq. (|D9[) have a factor of 2. Again, out 
of the eight possible combinations for s\, S2, and S3 only 
three combinations - Q( 2 ' ±+ \ Q( 2,± ~\ Q^ 2,t_ ^ - have a 
valid or non-zero analytical structure. 

The diagram Q^ 3 \ shown in Fig.[lOl describes the con- 
tribution of the correlation function with the interaction 
part of the heat current operator, J^' 1 ^ (see Appendix IB) , 
and has therefore a different structure in comparison with 
contributions Q'- 1 -' and Q^: 



g(3) = __L ei 2 T 



2^2 2^+2 ( (0) 



2d 



(Dll) 
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FIG. 10: (color online) Diagram describing correction to the 
thermoelectric coefficient due to electron-electron interaction 
corresponding to term in Eq. (|D2[) . All notations are the 
same as in Fig. [5] but here the heat vertex corresponds to j*' 1 
and the diagram has only one diffuson. 



with 

F (s lS2S3) = + ^ + q.jgtj.^ + n m )G?(oj n ), 

(D12) 

and 

53 (6, 6) = 2(l + ^t^) (D13) 

= 2(^i + ^-(a + 6))+o(e 2 /4)- 

Since the linear part of function (73(^1,^2) has a factor 
I 



Ep , the main contribution to from the diagram 
is of the order of T 2 /e F , whereas and have 

1/ef contributions, which we are considering here only. 
Therefore we will not consider diagram any further. 

In the following we discuss the five diagrams contribut- 
ing to r/ 1 ) in details, especially their analytical structure 
and the resulting restrictions on the frequency summa- 
tions: Q^'TT), Q(2>±+), Q( 2 >±~), and Q^~K 



1. Calculation of contribution Q (1±±) in Eq. (fD7l) . 

Here we discuss the contribution Q^ 1 ^^ introduced 
in Eq. (|D7|) . The analytical structure of this diagram 
(Fig. [5]), defined by indexes Si to S4 in Eq. (|D7[) . demands 

w„ + fi OT > , w„ + fi m + fij < (D14) 
w n > , w„ + fij < . 

These inequalities define the limits of frequency summa- 
tions: 

< lu h < — fi m — fii and fii < — fi m . (D15) 

First, we calculate the ^-integrals in Eq. (|D7|) using the 
residue theorem: 



/ 



a + a + f^ 



(*(w„ + On) - 6 + */(2r)) (z(w„ + fi m + fii) - 6 - »/(2t))(kj„ - 6 + V( 2 "r)) (i(u„ + fi 4 ) - £2 - */(2r)) 



-(2tt«) 2 
4tt 2 



z (2w„ + n m + (2uj n + fi m + fii) 



6,o=»(w„+n m +ni)-i/(2T)ec 
2 



4tt 2 



(fii - 1/r) 



(-tfi< + */r) (ifij - */r) (fij - 1/r) 5 

^ff (1 (2w„ + fi ro + fii)) « 4tt 2 t 2 5 (» (2w n + fi m + fii)) , 

I 



?2,o=iw„+i/(2r)eC+ 

i (2w„ + fl m + Qi) - - — (2w„ + fi m + fii) 2 



(D16) 



2. Calculation of contribution Q (1 ' T=F) in Eq. (TDT)) 

Here we discuss the contribution Q'-'-'TT) defind ir 



where in the last line we introduce the function 

g(z) = z + [d/(2e F )}z 2 . (D17) 

„, . , ,, c , ,1 Eq. (ID7I) . The analytical structure of this diagram do- 

lor the final approximation we used the tact that , — ' 

m £i ti ri ^ 

I fii I r <C 1. The poles for the residual are written be- 
low the under braces. 

Using the result of integration over a an( i £2 in 
Eq. (|D16[) . we can simplify the expression for this dia- 
gram to 



j n + fi m < 
uj„ < 



uj n + fi m + fii > (D19) 

LU n + fii > , 



i(i,±±) 



27T 2 S 



et 2 T 2 r 2 a 2d+2 (^ 0) ) 2 (D18) 



which defines the limits of the frequency summations: 

— fii < u> n < — fi TO and fii > fi m . (D20) 
We first perform the aintegrals m Eq. (|D7I) : 



0<^„<-n m -f2i 
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/ 



d£ide 2 (fi+e a +i^ 



i(u n + fi TO ) - & - i/(2t)] [«(w„ + n m + Qi) - & + z/(2t)][zw„ - & - */(2t)] [z(w„ + ^) - 6 + V(2r)] 

v . ^ . - 

£ 1)0 =j(u> n +Q m +fi i )-Hi/(2T)eC+ S2,o=«w„-i/(2t)GC- 
\2 



47T 



i (2w„ + n m + n*) - (2w„ + n m + fi<) 4n 2 g (i (2w„ + n m + 



(— — t/r) (iCli + i/t) 



(fi* + l/r) J 



4ttV<7 (» (2wn + n m + fi«)) (D21) 



Substituting the result of Eq. (|D21|) back into Eq. lfD7]l 
we obtain 

Q(^T) = _^f ei 2 T 2 T 2 a 2d+2 ^(0)^ (D22 ) 

x ^ ) (Cl i )g(i(2cj n + n m + Cl i )). 

— t>^n< — n m 

3. Calculation of contribution Q (2 ±+) in Eq. (fD9l) 

Here we discuss the contribution Q( 2 ' ±+ ) introduced 
in Eq. (|D9j) . The analytical structure of this diagram 
demands 

uj n + Q m > , uj n + Q m + Qi < 

w B > (D23) 

which defines the limits of frequency summations: 

< w„ < — fi m - fij and < — fi m . (D24) 

I 



For symmetry reasons we write both versions of the in- 
tegrals (diffusons on grain 1 & 2) and therefore already 
take the factor 2 in front of the subdiagrams into ac- 
count in the right hand side of Eq. (|D9[) . Furthermore we 
introduce the short notation for the Green's functions: 



o±(0 = i{u n +n m )-t±i/{2 T ), (D25) 

c±(0 = iw n -f±*/(2r). 



In the following we need only the poles of functions b±(£) 

and c±(£): ^ = i(w„ + fi m + fi f ) ± i/(2r) G C ± and 

= iu„ ± i/{2t) € C ± . We also use the function 1712 
introduced in Eq. (|D9[) . 

Therefore the ^-integrals in Eq. (|D9j) can be written 
as: 



/ 
/ 



^&<?(a+&)[(4(^)M6:M&)) 1 + (4(&)&_(k)c+(6)) 1 

«ide 2fl (6 + 6) a2+(6 )6_(6)c + (a)a^(6)6-(6)c + (6) 



2m / c^i 



<(£i)M£i)c+(6) 



2m / d£i 



2<(a)6-(a) 



4 (d b) 



(d b >+(a) 



(D26) 



Here we executed the ^-integral and used the fact that 
the first term in the second line has only one pole in C + 
(factor ni). We are left with two term where the second 
one has only a single £i-pole in C + which gives another 



factor m. However, both terms give the same result 

~(d c) +d b)N 



4tt^ 



4tt 2 5 (2w„ + n m + Qi)) 



-4ir 2 T 2 g {1 (2w„ + Q m + f^)) 



(D27) 
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Substituting the result of Eq. (|D27p into Eq. (|D9j) we 
obtain 

2Q (2,±+) = ei 2 r 2 r 2 a 2d + 2 ^(0)^ 2 ^ 

x ^ ^ (n l )g(i(2co n + n m + n l )). 

o<oj„<— n ra — f2j 

f2 j <— n TO 

The notation introduced in this Appendix allows us to 
write down the Q^ 2 ' T ~^ and Q( 2 ' ±_ ) contributions in 
Eq. (|D9p by just changing the +/- indices. 



demands 



uj n + O m < , w„ + fi m + Oj > 0, (D29) 

Un < 0, 



which defines the limits of frequency summations: 



4. Calculation of contribution Q (2 T_) in Eq. (|D9|) 



O m — f2j < cl>„ < — TO and f2j > 0. (D30) 



Here we discuss the contribution Q^ 2,T > introduced 
in Eq. (|D9|) . The analytical structure of this diagram 



Using notations introduced in Eq. (|D25[) the ^-integrals 
in Eq. (|D9[) can be calculated as: 



^i^tfi +6) [(aL(6)M6)c-(6)) 1 + (oi(6)6+(6)c_(6))" 
, ^ a 2 _ (6)6+ (6 )c- (gi ) + a 2 _ (gi )6+ (6 )c- (6 ) 



(D31) 



2m / d£ 



a 2 _(Ci)& + (a)c-(a) 



.9 




a 2 - (6)6+ (fi) 


a 2 _ 







2m / dfx 



2« 2 -(a)&+(a) c_&) 



4tt z 



(d b) ) 



-4tt 2 t 2 3(i (2w ri + n m + Oi)). 



r 



Substituting the result of Eq. (|D31|) into Eq. (|D9[) we demands 
obtain 



2Q (2.T-) = ^ et 2 T 2 T 2 a 2rf+2 ^( ))( D3 2) 

x £ *£ 2 >(n i ) 5 (*(2w„ + n m + n i )). 

n i>0 



5. Calculation of contribution Q* 2 '*"' in Eq. (|D9[) 



w n + n m > , u n + fl m + Qi < 0, (D33) 
ui n < 0, 

which defines the limits of frequency summations: 

- il m < uj n < and fi, < -fi m , (D34) 
and the disjunct region 
— fl m < Lu n < — Q m — Cli and — fi m < f2,- < 0. (D35) 



Here we discuss the contribution Q^ 2,± ) introduced 
in Eq. (|D9|) . The analytical structure of this diagram 



Using notations introduced in Eq. (|D25[) the ^-integrals 
in Eq. (|D9p for this diagram can be calculated as: 
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JdZtdZtgfa + b) [(a^(a)6-(a)c-(6)) 1 + K(&)6_(f 2 )c_(^i)) 1 

MM u , , v 4 (6 )b- (6 )c_ (6 ) + 4 (6 (6 )c- (6 ) 
d£,id£,29 (£i + 6) „ 2 ^ TZ ^ ~ ^ \„2 — 



(D36) 
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2tti / d£i 
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2a 2 



(6) 



(a)6-(6) 4 (d b) ) c _(ex) 



-4tt 2 



,. 9 (d c) + d b) 



4 (d b) ) 



47T 2 T 2 ff (t (2cj„ + Q m + fii - 1/t)) 



Substituting the result of Eq. (|D36[) into Eq. (|D9|) taking into account the two disjunct regions for the frequency 
summations, we obtain 



2Q (2,±-) = ^l et 2 T 2 T 2 a 2 d+ 2 ( (0)- '< 



E + E 

-fi m <u> n <0 -O m <o; n <-n m — . 



(fii) .9 (» (2w n + Sl m + (U - 1/t)) . 

(D37) 



6. Analytical continuation 



Here we combine all five contributions to calculate correction 77^^: (5^ 1,±± - ) , Q' 1 '^' , Q^ 2,±+ \ Q^ 2 ' T ~\ and Q^ 2,± ~^ 
introduced in the right hand sides of Eqs. (|D7p . (|D9|I and discussed in the previous subsections (|D ip through (|D 5p . 
We focus in particular on the analytical continuation of the Matsubara to real frequencies, both for the Fermionic 
and Bosonic frequencies. Using Eqs. (|Dl8]) . (|D22j) . (|D28|) . (|D32]) . and (|D37f one can write 



QW + Q( 2 )=Q( 1 -±±)+Q(^T) 



2 Q( 



(2,±+) j_q(2,T-) j_n( 2 .±- 



AT 



T 2 ^(n i )g(i(2u )n + n m + n i ))- ^ r 2 $« (O i ) 5 («(2a; n + O ro + O i )) 



o<u;„<-n m -ni 

Sl i <-Sl m 



-Qi<u)„ <— n„ 



t 2 ^ {Q t ) g (i(2uj n + Q m + Q t ))+ t 2 ^ {a i )g(i(2uj n + n m + Q t )) 



0<UJ n < — Qm—Qi 
fJ»< — (lm 



-f2 m -Qi<a; 7Z <-n, ) 
fij >0 



E 

-n m <w„<o 



r 2 $<^ (fi<) ff (2w n + O ro + n { - 1/t)) 



E 

-n TO <w n <— n m — 

-S"2 m <f2i <0 



(D38) 



- 2 $w (a) g {1 (2w„ + n m + fij - 1/r)) 



where we introduced the notation A = 2jI j^et 2 l Ta 2d+2 \^d) anc " tnc functions (Hi) and <f>£ 2 ^ (0^) were defined 
in Eqs. (|D4|) and (|E>5|) respectively. 

Next, we perform the summation over Fermionic frequencies u n in Eq. (|D38|) by shifts and the corresponding 



analytical continuation 
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o<a;„<-n m -ni 

n,<-s7„ 



^2 /(w„,0. 



-ni<w„<-f2„ 



-f2 m — fii<(D„<-n„ 

n, >o 



f(w n ,Qi 



-n m <w„<o 

!2;<-n m 



-Q m <ui n <-fl m -Qi 



E E 

w„>o u„>-n m — Oj 



/(w„,a)= V [/K,0^)-/(a; T[ -f2 m -f2 i ,f2i)],(D39) 



Ci i <-Ci 7 



E E 

a; n <— Q, m u) n < — Qi 
n i >n m n i >n m 



SI,- >£i m 



E E 



/(a; n , n*) = y2 [/(w n - n m , n<) - /(a; n - - n m , n^] 



E E 

sii<-fi m n,-<-n m 



/(w n ,f2t)= ^ [/(a; n ,Oj) - /(a;„ - , 



E - E 

;„<— fi m — f2j u) n <-Q 

= [f(w n -ft m — Sli,Sli)-f(w n -Q m ,Qi)] 

Ul n <0 

—Q m <Q i <0 

= J~] [f(u n - Ct m - Clj, Slj) - f(w n - Ct m , Ctj)] 

UJ n <0 
f!;<0 

- ^2 [f(u n - TO - Cli, fi 4 ) - f(uj n - fi m , fij)] . 



OJ„<0 
S2 i <-S2 r „ 



Here the function f(oj n ,fli) is the product of the func- 
tions t 2 $L q ' ) (^i) and g. For the analytic continuation we 
need to consider the to -dependence of function $>ffl (Oj), 
in particular we can use the fact that (which follows di- 
rectly from the definition of functions in Eqs. (|D4|) 
and (jD5|» 

^ = ^i (D40) 



where $g is the same as $L but with r w replaced by 
To and is therefore w-independent. Using Eq. (|D40|) , we 
obtain 
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XT 
~ AttiT 



Q(D + Q(2) = Q(l,±±) + Qd-TT) + 2 (QP.±+) + Q(2.=F-) + Q( 2 .±-)) 



(D41) 



(^)-4 2) (^) 



rfw tanh (^;J [5 + * (fi m + fij)) - g {2lj - 1 (Q m + CI,))] 



+ J] ^ (Oi) To / rfw tanh (^) [9 (2w + t (-fi m + n 4 )) - .9 (2w + * (O m - Oi))] 
- ^ $< 2) (Qi) r 2 / ^ tanh (|M [g {2u + 1 (-Q ro + Q;)) - .9 (2w - 1 (O m + SI;))] 

f2»>0 



£ $( 2) (Qi) r 2 f dcu tanh (^) [.9 (2w - 



! TO + 0i-l/r))-ff(2w + «(- 
! TO + Oj + 1/r)) - 5 (2w + 1 (-0, 



^77 



i m + n< - i/t))] 

n,<-Q m J 

- $ o 2) (^) T o I ^tanh (|M [g (2uj - t (ft m + SI, + 1/r)) - ff (2w + 1 (-f2 m + ^ - 1/r))] 



2 4 2) (0*) r 2 / du tanh [g (2w - 1 (fi m + 0* 



- 1/r)) - g (2w + t (-fi m + n t - 1/r))] 



Next, we consider the integrands in Eq. (|D41|) and introduce the short hand notations 



01 = .9 (2a; + 1 (0 TO + fii)) - 5 (2w - i (O m + Vti)) 
a 2 = .9 (2lu + 1 (— fi m + fij)) - .9 (2w + 1 (O m - fi;)) 
a 3 ee g (2w + ? (— fi m + Oj)) - .9 (2w - i (O m + fi;)) 
a 4 = .9 (2w + 1 (f2 m + Oj - 1/r)) - .9 (2w + 1 (-fi m + ft, - 1/r)) 
a 5 ee .9 (2u - 1 (f2 m + + 1/r)) - .9 (2u + 1 (-fi m + £1, - 1/r)) 



2z (f2 m + 0i) (l + 2dw/£ F ), 
2i(fi* - fi m ) (1 + 2duj/e F ) , 
2z^i + 2£l;d (2oji + tt m ) /e F , 
2Q m t + 2n m d/ep (1/t u + 2ilu — fij) , 
-2iD,i - 2dQ.il e F (l/r u + 2iuj - fi m ) 



(D42) 



r 



Now, we extract only terms which are linear in lu in 
Eq. (|D42[) and of order 1/ep (thus the r w in 04 and 05 
does not give a contribution). Therefore we obtain 



ai ~ 2lji (Q m + f2,-) 2d/ep 

a,2 — 2lji (fi, — f2 m ) 2d/s F 

as ~ 2ujiQ.i2d/sF 

~ 2unSX m 2djEF 

a5 ~ —2un£li2d/£F- 



(D43) 



Substituting the result of Eqs. (|D43)) back into Eq. ([D4Tj) 
we obtain 
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2r 2 z2dAT 



Qi<-a 



(1) + g (2) = q (1,±±) + g (1 - TT) + 2 (y 2 < ±+ ) + g (2 ^-) + q(s.±-) 

- ]T (0^) - ${, 2) (fi<) A tanh(a;/2T)w (Q TO + Q;) 

n;<-n m ^ 

^ $^ (fij) / dwtanh(w/2T)w(f2i-n m )- ^ $£ 2) (fii) / dwtanh(w/2T)wO i 

f2i>n m ^ o ; >o 

^ *o 2) (^) / dwtanh(w/2T)wO m - ^ $[, 2) / dwtanh(w/2T)wO i 
^ <4 2) (n 4 ) A dwtanh(w/2T)wf2 i 



fli<0 



Qi<-n 



(D44) 



[ $ 1} - $ 2) («<)] + Hi) + J! $ 1} - M - E $ " 2) ^ 



f2»<-n„ 



a<>n„ 



^ $< 2) a„ - ^ ${, 2) (no a + ^ $< 2) (fij) a 



f2j<-a„ 



f)i <o 



n t <-n„ 



In the last two lines of Eq. (|D44j) we introduced the no- 
tation A = J dwcJtanh(w/2T) where we use the fact 

that J duju tanh(cj/2T) — >= — , neglecting the in- 
finite boundary terms of the partial integration. 

Using the fact that the functions <J>^ 2) (fij) in Eq. ([D44]) 
depend only on the absolute value of the Bosonic fre- 
quency Q, we finally obtain 



Q«+Q< 2 
= -2Ai 



E 

n»<-n r 



(D45) 

$^ 2) (n m + n«) , 



where we 



use 

,(0) n d+l 



the 



notation Ai = 
The internal frequency 
summation is also done by analytical continuation, but 
for Bosonic frequencies: — > —iQ + r\. Before the 
final integration we should take the external frequency 
derivative and finally calculate the q-integrals of the 
<l>o-functions. Using Eqs. (|D4[) and (|D5[) we have 



(2) 



dq 



2E c (q) 2d - J2 a cos (q ' a ) 



I" 2 (|^|+ 4^(9)^) (|^ 



Using Eq. (|D46[) one can calculate the sum over internal 
bosonic frequencies f2j in Eq. (|D45() 



Sn = 



E 



(D47) 



Qi<-Q, 



E 



f2i<-n, 



E 

Oi<0 



(|nj|+4E c (g)e,)(|ni| + e,<y) 

(-^+4S c (g)e g )(-a i + e g 5) 

Sk 

(-Qi + fi m + 4£ c (<j)e,) (-fij + n TO + e g 5) 

[-4 7 rzT]- 1 (-zf2) coth(f2/2T)dft 
[*fi + ft m + 4£ c (g)e g ] [iQ + Q m + e q S] 

[-47rT]- 1 f2coth(f2/2T)dfi 
[O - t (n m + 4B c (?)e,)] [A - t (O m + 

[-47rT]" 1 a;coth(a;)dx 
[a: - ^ (fi m + 4E c (g)e,)][a: - ^ (fi m + e q S)} 



(D46) Only the terms proportional to the external frequency 
TO in Eq. (|D47[) contribute to the correction to the ther- 
moelectric coefficient r/ 1 ) in Eq. (|D2[) . Taking the deriva- 
tive of both sides of Eq. (|D47p we obtain 



2G 



d_ 

on 



Sr 



1 



4irT 
1 

~8ttT 2 
1 



d_ 

on 



x coth(x)dx 



n 

2T 



a; coth(x)<ix 



1 



n 

2T 



(D48) 



(x — ia) 2 (x — ib) (x — ia) (x — ib) 2 



x coth(a;)c?a; [2x — i (a + b)] 
(x — ia) 2 (x — ib) 2 



8ttT 2 



Ia,b, 



where a — 2E c (q)e q /T and b = e q S/(2T). Finally, we use 
the approximation: 



x coth(x) 



1, |ar|<l 
\x\, \x\ > 1, 



(D49) 



and obtain for the integral I a ,b in Eq. (|D48[> the following 
result 



I a .b = i 
For a > 1 > b 



ln[(l + q 2 )/(l + fr 2 )] 
(a-b) 



2 In a T\n[2E c (q)e q /T] 

la.b = i = t jrn • 

a E c {q)e q 

Thus, for Eq. (|D48|) we obtain 

-* ]n[2E c (q)e q /T] 



d_ 



Sq 



n=o 



8ttT E c {q)e q 



(D50) 



(D51) 



(D52) 



Substituting Eq. (|D52|) into Eq. (|D46j) we obtain for the 
correction to the thermoelectric coefficient the following 
result 



A 



x / dq- 



2nT 2 (2TT) d T$ 

2d - £)' a cos(q ■ a)l In [2£ c ( g )e q /T] 



(D53) 



7. Final integration and results 



In Eq. (|D53[) we are left with the final integration over 
internal momenta q. Therefore we need the functional 
dependence of e q and E c (q) on q, which were introduced 
around Eq. (|D6|) as 



The final q-intcgral can therefore be written as 
-A 



V 



(i) 



, dqln[2£ c (<z)e q /r] . (D56) 

47TJ Z (2TT) a gTTQ ' 



The q- integral is cutoff at q = n /a. In d = 3 the In - 
argument is finite at q = 0. In d = 2 the volume element 
makes the integrand at q = finite. The coefficient in 
Eq. (|D56j) can be simplified to 



-A 



ir 2 s et 2 a 2d + 2 



47rT 2 (27r) d ffT r 2 6(27r) d g T 
fj{0) / a \ d 



VP) 



2 T 

e F 



such that in d - 
,(<>) 



2ngT V27T/ 
3 we obtain: 

3 



(D57) 



AY 



x In 



2irgT \2nJ 
lQite 2 gT 3 — cos(aq x ) — cos(aq y ) — cos(aq 



(D58) 



aT 



'/ 



(0) 



2ng T 



f a \ 3 7r 3 
V2W ^3 



M 2 

4tt, / e 2 

y m I 16— gT-n/T 



C3 



Here C3 = / dqln (3 — cos^^) — cos(7rq , j / ) — cos(7rg z )) / (irq)' 

and integration is over the unit sphere. However, this 
numerical constant (and all other inside the logarithm) 
can be neglected, since E c /T S> 1 and gr ^> 1. 
For d = 2 we obtain 



Co) 



2irgT V 27r 
8e 2 <7T"" 



with C2 



x In 

_^(0) 
27T.9T 

= f ^ In 



To 



rfq (D59) 
(3 — cos(ag K ) — cos(agj,)) / (og) 



/ a \ 2 7T 3 
V2W ^2 



In 



Ta 



C2 



2 — cosfTrqx) — cos(7rq y 

77^ 



2<? T 2d - J2 cos(q ■ a) 



(D54) 



E c {q) 



2C(q) 



— ]n(qa), d = 1 

ir/q,d = 2 . (D55) 

2Tr/q 2 ,d = 3 



a. Final results 

Now, we can summarize our final results in d = 2, 3 
for the thermoelectric coefficient (only the functional de- 
pendencies of the correction terms are kept under the 
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logarithms): or combined in a compact way (valid for d — 2,3) we 

obtain the following final result for the thermoelectric 



•S7T 



3 



2 d+2 ( (o)\ 2 T (r,an\ coefficient of granular metals 



„<°> = ^—et^^>) f-, (D60) 

(i) S7T 3 ea 5 2 / (0 )\ 2 T ( e 2 g T \ m . 

= m (e c9t/ t) , „ = ,(o) (i - ^ m ^ ) . (D63) 

- ^-M 0) ) 2 -mf^ (D62) 



2A g T \ J £f \Ta 

„(0) 

■i- \n{E c9T /T) 
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